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The ccmguter program calculates tne supersonic flcw over 
a cere in air or water. The rain objective is to calculate 
the ccne semi-vertex angle given prescribed initial ccndi- 
tions. The program is written in structured FCRTEKAN and 
itplezents Busegzann's graphical integraticn technique. 
Superscnic flow over a cone in water is useful as a gcod 
first appreximation tc the motion of the metal jet frem an 
explosive shaped-charce fired underwater. 

A typical result for supersonic flow over a ccne in 
water is as follcws: given an upstream temperature, 323.16 
Kelvin; upstream pressure, 1 Lar; shock angle, 20.0 degrees; 
and pressure behind tke shock front, 5 kilokars, the ccne 
semi-vertex angle is calculated to be 7.23 degrees. 

Generally, Eresstres invclved in water flow are much 
darger tkan for air flow, and the cone seai-vertex angles 
for water flow are smaller than for air flcw. 
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I. INTEOLUCTION 


aA. CESCEIEITION CF TEE FROBLEM 


The sclution of tae hydrodynamic equations describing 
superscric flow cver a cone in air has teen well kncwn since 
the 1630's. Until recently, the problem of descritince the 
flow cver a cone in water has been limited to scluticns of 
the sukscric case. Frimarily, calculations were lisited to 
subscric flcw Lecause researchers believed that supersonic 
flow in water was not feasikle for normal vessels (such as a 
ship). Fer ordinary vessels in water, it is certairly true 
that superscnic flow past that vessel is highly imgprckable. 
HOWwEVEL, the mction of the metal jet from an explesive 
shaped-ckarcse fired ctrderwater is supersonic. 

This thesis ;reserts a ccmuputer program whicn calculates 
the hydrcdyramic flcw fast a cone ia either water ocr air 
under stpersonic ccnditions. The program utilizes the 
methcds developed by previcus researchers fcr calculating 
the superscnic flow in air and which have been = suitatly 
modified tc describe the conditions in the water. Suca 
nodificaticns include utilizing the modizied Tait eguatica, 
which is the "thermal" cr "thermodynamic" eyuation cf state 
for water, to describe the j;hysical state of the sater 
rather tkar the serfect gas law used for air. 

In actuality, tke cone lirer in the jet from an explo- 
sive shaped-charge is blunt-nosed rather than ccrical. 
However, solution cf the ccnical case is a preliminary 
requirement to soluticn of the actual blunt-nosed fretlen. 
The scluticn to the conical flow case will serve as an 
excellent test prograzr for the solution to the blurt-rosed 
Frobler. This thesis presents a solution to the ccrical 


case anc it is aojed that the projran wiii assist the 
contiruirg research into tne froblems of utilizing exylcsive 
shaped-clarces in ar tnderwater environment. 


Ee. HBETECDCICGY 


Ike ccmputer prcgram presented in this work was origi- 
nally developed in the BASIC computer language usirg a 
Hewlett-Fackard HP-67 computer. That program is the kLasis 
for this thesis. It was desired to translate the fregran 
into a tEigker-order computer language for executicn cn a 
larse, wainframe cciputer systen. This translaticr vxas 
desired in crde. to gake the program more readily accessitle 
to a wider pody of researchers and in order to sfeec the 
executicr time of the prcgran. In this thesis, the 
fecllcwing gcals have teen accogrplished: 


e Successfully trarslate the program from BASIC into a 
higker-crder lancuage. This joai was met by utilizing 
FCEFIFAN as the tign-level languaye of choice. khile 
FCBHIRAN has gany drawbacks as a high-order language, it 
is still widely used in the scientific ccmmunity. 
FCETIFAN was used, therefore, so that the frogram will te 
useful to as wide an aré€a cf researchers as pessitfle. 


e Fellcw modern picgramming g¢ractices in the design and 
iaoplemertaticn of tne program. As before, the choice of 
FCRIFAN as tke high-level languaye makes this gcal some- 
what wore difficult. However, Many computer scientists 
have demonstrated that structured proyramming practices 
can té achieved tusing FORTRAN. To the largest extent 
Fessifle structtred proyramminy fractices have teen 
utilized. 


e Fresent a "user-friendly", well-documented prograr. In 


this regard, literal use of comments occur in the 


CPM GF ates mo 


ELroglam itseli, geaningrul variable names are used, and 
detailecé flowcnarts whici demonstrate tae logic cf the 
Eregrén are incltded. In addition, due tc iitited 
interaction witk the user, the user's respenses are 
verified befcre tke prograg executes. 


1C 


Similarly, ir cortrast to air, the thermodynamic Changes 


woick cccur as a result of the saock process in water canrot 
ke e€asily delineated Ly simple equations as in the air case. 
However, a simplification can ce made in tne water case 
kecause, urlixe in tke air case, tne pressure jump acrecss 
the skcck in water is so very large. In air, [resstre 
changes acress the skcck f£rent cn the order of 1 toc 2 bars 
are ccnsidered large {at least for chemical explosicrs). In 
contLast, as pointed cut by Richardson, et. ale, [Bef. 3], 
the fFressure jump acrcss a shceck in water is cn the crder cf 
kilotars te tens of kilobars. Therefore, the calculaticns 
can Le simplified ry specifyirg the pressure on the écwn-~ 
stream side of the skcck frent. This is valid since the 
upstream pressure is so small in comparison to the upstrean 
dynakic fressure Oo; vi 72 and ain comparison to Ecth the 
dcewnstream pressure ard dynagzic pressure. The specificaticn 
cx the dcwrstreag pressure is accomplished, in the precgraa 
cr this thesis, by allowing the user to ingut a "f;resstre 
gultiplicaticon factor (MFACI)", which converts the j;resstre 
upstream tc a pressure downstream at pcint 2 which is given 
ky: 


E, = F, X MFACT ¥ 1000.0 (2.24) 
Toe tactcr 1000.0 ir eguaticn 2.24 converts the right-hand 
side cf tne eyuaticr from a pressure in bars to a [fresstre 
in kilckars. As an example, if F, is 1 bar and MFACT is 5.0, 
the pressure at foint 2 downstream will be 5 kilobars. 

It can be shown that the Simjiifyiny assumpticn made 


above is ertirely valid by considering the momentum e€guaticn 


for steacy frictionless flow alcny a streamline. The 


owe 


Te. ECCATICAS SPECIFIC IO HATEE 


Ir ccnotrast to air, which has a relatively elegant a ni 
Simple state e uation, the eguarion of state for water is 
Latner acre complicated. The ncst commonly used state equa- 
tion for water is kncwn as the modified Tait eguatior, which 
may Le written as: 


o=BS)( - 3] (2-21) 
p 


shere F(£) is a slowly varying function of entropy alcne, n 
is apprcximately a ccrstant equal to 7.15, and ¢ is the 
value cf the density for zero pressure. The above eyuaticn 
is frcem Eclt [Ref. 6}, but, in different forms, it is also - 
descrilked guite extensively in Cole [Ref. 2], kicharéscn, 
€t.._a1., [bef. 3] ard Kowlinscn [Ref. 7]. AS merticred in 
the Irtrcduction to this work, the modified Tait eyuaticn is eae 
the "thermal" egquaticn of state for water. Cole [Ref. 2] _ 


shows that the modified Tait equation is of the form 


n as 
: = 3s) (ED eT, (<. 22) Doe 
cr, ip sigpier terms, 
F = plv,T) (Z. 23) 


Since tke full modified Tait eguation relates the three 


therucdyragic guantities of ¢~, v, and T, it is calleé the 


"theracdynansic" cr "thermal" equation or state. 
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Figure 2.4 Gecasetry of the Obligue Shock Front. 


2 eed k ~ 1 
Pp kM, sin’B Coa) 
1 (——) 


Eqguaticn 2.17 is used to determine the pressure dcwnstrean 
cf the stock frort. 
k - 1 2 0.2 / re k - 1 
T, c, 2 1 + ( 5 )(M) sin E) (kM) sin’ - [ a }) 
ro eo a oe ee I Oo eo cage pe ee ae ee 
1 1 ( My 


k + 1 
2 


) (Z. 18) 


sin 6 


Equaticn 2.18 is useé to determine the emperature, and acre 


importantly, the speed cf scurd, c , downstream of the shcck 


front. 
u u,sinf * 
M =—s pe ee = My sing (Zz. 19) 
x Cy cy 
u u,sin(B - 4) AS a 
M = —2 = 2 = M,sin(B - 6) (2.20) 
y Cy Co Z 


Finally, é€guaticns 2.19 and 2.20 are used tc determine the 
velocity componerts cf the flow across the shock frcrt. 

The cecnetry of the flcw conditions across the shcck 
front is illustrated ty figure 2.4. Note that in the e€qua- 
tions akcve and in ficgure 2.4, the *1* subscripts refer to 
conditicrs in the trshocked (upstrean) fluid, the ‘x! 
subscripts refer to the normal components cf flow in the 
unshcecked fluid, the '2' subscripts refer to conditiors in 
the skccked (downstream) fluid and the ‘y' subscrifts refer 
to the nergal components of flcw in the shocked fluid. Ncte 


also tkat *k* in the eguaticns above is the designaticr for 
the ratic cf the heat capacities Ci/Cy° 


a 
te 
a 
a 
aa 


wosre Ris tae specific sas Hetant and as reiated te the 


ce 
universal gas corstart, \, cy: 7s 


» 
a 
<|> 
x 


(<- 14) 


In egtation 2.14, M, is the gclecular weight of the air. 

Ir aiz, the charge in the thermodynazic progerties of 
the ga¢ as it cresses the shcck front are easily calculated, 
as Shcwn io Kinney ard Grahag [Bef. 5]. Since most fluii 
dynasics teatbooks illustrate the develo, ment of the ¢€qua- 
tions abich folicw, it is net necessary to derive them here. 
As werticned previously, Kinney and Graham [Ref. 5] grevide 
éxcerticnally lucid explanations and derivations. MTke frin- 
ciple eguations used to calculate the thermodynamic ckances 
thich cccar across tke shock front in air are as follcws: 


2,2 
tan(2 - 6) | 2+ (k ~ 1) M, sin 8B (z015) 
ne (k + 1) Mi sin?e soo 


Equaticn <.15 is used to iteratively determine the deflec- 
tion angle 6. All cther guantities in this equaticn are 
known (i.e. 8 is tke snock aajle and a is the freestrearn 
Mach ouster, beth cf which are input ;arameters to the 
Erogras for the air calculaticns). 


> 2ttk-1 Mi tard 
[M,sin(B - @)]° = —3>——;———— (Z. 16) 
2k sin’B - (k - 1) 


baving determined 6 from eguation 2.15, equation 72.16 is ide 
tsed tc determine the Mach nunter on the dcwnstream side of : 
the sktcch front. 


- ~~ diy, MEE KBEN 4 Sua eho hS . 
i rr ce a Rote te we ote ee ee a te ee et 
ete te te fe te oe 8 te hg Me Me Me ee Oe Me Og ee es OP aes 
Po ee Pe el a al Se a 


oat te tea te te te oe ot es eh ee eh te Mee et et 


Py ne ae ne ° ae Lt ta ee Tle) 


Oe a ne Pa le Lee ete ee eee ae ee as 


Ey ccacirirg Zqn 2.9 with Egn z.8, one arrives at: 


sin = 
V_ sinw (2-10) 
a eae: | 
2V" sin (w - 9) 
ae 3 3 
(k - BY Vs ah¥-s) 
Eut, the energy eguation asserts that: 
2 k - 1 2 2 - 
= = - 11 
c es Pn ED) (< ) 
therefore 
sind 
R = V sinw 
ve sin? ( - 8) (2. 12Z) 
l-——_—_ 
2 
Cc 


where c is the lccal speed of sound in the fluid. 

Eguaticn 2.1Z is the basis for the calculaticn cir the 
superscric flow cver the solid cone in either air cr water. 
The c¢xaphical iategration method invented by Buserann is 
adequately explained in Shapiro [Ref. 8] and need nct be 
repeated Eere. Essertially, the computer proyram given in 
this werk automates the Eusemann graphical integratior 
metoce fcr calculatinc the ccne semi~vertex angle. 


Cw kECUATICNS SPECIFIC TC AIR 


The e€gtaticn of state fcr air is specified Ly the 
Ferfect gas law (under the assumption, that is, that the air 
behaves as a perfect gas). This law is quite elegant and 
Simple ard allows easy manijulation to obtain various guan- 
tities. The form of the perfect gas law most often used in 
the calculations of this thesis is: 


p = oft (Z. 13) 


s EM ee eat Eta S Wed Nid Pah tay Feel, ne am Ui Bh, ee aN eed ie ici ea Vaal aD 
Teg) cet Te oe Ser th pie eg eee artic Sp ee ee Mee we r 
ea my? ete Me” 7 = = en Te Me Sw et ’ x 


a 


Eut, fren Egn 2.1 and Fgn 2.2, it can pe seen that: 


=V = -V sin(w - 6) (2.5) 


_ WwW (2.6) 


dw : Vtan(w - §) 


Shapiro [Ref. 8] shows that the eguation yoverniny the flow 
in the ccrical regior is: 


dv 


k - 1 Ww 2 ee 2 
( 5 )(2V_. + vy cotw + mR VS - ve - Vv) 
(2.7) 
dv 
r dv 
= We dae Ms sv 
iw’ w 


Elinirating * aber Egn 2.4 and substituting the exyres- 
; Vv dv . a 
sions fcr ve Vo P = and at Given by Eqns 21, z.4 and 


z-5, inte Egn 2.7 gives: 


sin® 
dv Vsin(w - 8) sinw 
dea 7 ON statu - 9) (258) 


l 2 


2 
OS EN ee ON) 


Tesignating ‘k* as tke radius of curvature of the hcdccraph 


streanlire, one cbtairs: 


pe Seeeets:) Seareaee (2.9) 
sin(w - 8)dw 


18 4 
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In tre iolicwing derivaticnos, as per Shapico [fef. €], 
tae sjrericéel ccordirates r and . aave been used witr the 
ccrres;crdirng velocity comgernents we and Ne (see figere 
z.3). Coly tne prinary eguations wnico are used in the 
computer frcgram are presented in tais thesis. A détailed 
derivaticn cf the e€quaticns can be found in Shapiro 
(Ref. &j, and need nct be repeated here. The nomenclature 
used in the equations develcsed in this chayter is detailed 
in tTatkle I. In keeziny with modern thought, the M-k-S 
(metéer-~kilcgram-secord) unit system has been used thrcughcut 
this thesis except fcr occasicnal lapses during the water 
calculaticrs when pressures are referred to in urits of 
kilctars. 

Free the geometry of figure 2.3, it can be seen that: 


Vv. = V cos(w- 6) and an = -V sin(w - 8) (z.1) 


In the development ct the actual second-crder differertiai 
€guaticn, Stapire [Ref. 8] shcws that, due to the ccndition 
cf irrotaticnality, tke following relation must be true: 


dv 
Vv = pes (<.2) 


Differentiating E,n. 2.1 with respect to » , cne oktairs: 


dv 


ae GON tints av 7 (2.3) 
ae -V(l - qo sin(w 8) + do cos(w 8) 
and 
dv 

De = AB bya Se 2.4 
aa -V(1 FRG cos(w 8) = sin(w - @) ( ) 
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Figure 2.1 Shcck Cone and Typical Streamline. 
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After caiculatine the thermodynamic cnanyes which cccur 
as a result of tke skccx frcot, it is assumed that the fluid 
properties siil remain constant on imayinary cones kaving a 
commer vertex. By this assumftion, the flow past tne ccne 
can te calculated. The flow geometry is illustrated by 
figure z.1 and figure 2.2, which Show a typical streamzline 
and its izage in the Lodograph flane. 

As discussed thorcughly by Shapiro [Ref. 8], there is a 
disccrtinucus chanje in both direction and velocity acrcss 
the ckligue shock frent. Pcints 1 and 2 (see figure z.1 and 
figure 2.2) lie, tterefore, cn a common Shcck pelar which 
cliginates at point 1. Between point 2 and the point '"s', 
which is cn the core surface only at infinite distance, 
there is a regicn of conical flow where the stream fircper- 
ties vary ccntinuously. The velocity vector to the foint 
"st ain tke hodograpk flane defines what is called the ccne 
semi-vertex angle witk the centerline axis. In the methceds 
which fcllicw, the ccne semi-vertex angle is the variatle 
which is wltimately cetermined. Further, as pointed cut by 
Shapire [Ref. 8], since all streamlines in the flew exyeri- 
ence tke sane entropy jump acress the shock front, tke flow 
ketween tre shock ficnt and the cone surface is both isen- 
tropic ard irrotatioral. 

Ike second-crder differential eyuation which ectually 
descrikes tke flow cf the fluid past the cone in air is 
fully ceveicped wy Sktayiro [Ref. 8]. Shapiro aotes tkat 
there are two methods commonly used to solve this eguaticnh. 
Cne methkcd, developed first ty Taylor and Maccoll [Ref. €], 
perfcxrms a cumerical integratica of the equation. The seccnd 
nethcd, which intecrates the eguation using a grarhical 
constructicn method, was first developed by Busemarn. the 
Frogran developed in this thesis utilizes Busemann's methcd, 
ucdified fcr performarce on a rodern high-speed computer, to 
gerform the integraticn of the full second-order differen- 
tial é€quaticn. 
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ike funcameotal e€cuation used to describe tne teeracdy- 

. ramic state of air is the jferfect jas law. In water, the 

7 modified Tait eguaticn is the eguation most often used to 

éescrike the thermodyramiv state of the water. The mcdified 

F Tait €guaticn can be used tc describe either pure water or 

: sea water. The forr of the modified Tait eguation used in ° 

this thesis was taker frem Hclt [Ref. 6], whe has ccotinued 

to perfcrm research ir underwater explosion phenomena. Au 

excellent discussion cf the mcdified Tait eguation ard how 
it can ke utilized is contained in Rowlinson [Ref. 7}. 

After tre thermocynamic ;rcperties cf the water (cr air 
as tke case may be) on the dcwnstream side of the shcck 
front have feen calculated, an iteration method, utilizing 
an autcmated graphical ccnstruction first develcped by 
Eusemann, [kef. 8], is used tc progress from the shcck frent = 
to the ccne surface. The equations needed for use Ly this 
iteraticr method are fully described, for air, ky Skapiro 
{Ref. 6], who descrites their development and use. The 
Methocs which apgly tc the flaw rast a cone in air can, wita - 
the necessary changes made for the differences in tke ther- iS 
wodynarics of the two fluids, be used to calculate the 
conical supersonic flow in water. These caiculaticns icra 
the kasis for the sain part cf the FORTRAN pregrag shich ™ 
follcas. It shculd ke noted here that Shapiro alsc fgeinrts 
: cut the ;icreerinyg werk of Taylor and Maccoll in the 1630's 
7 and 1$40"'s on the mé€thods cf sclution of the flow over a 
a cone in air problem [Fef. 8]. 


r E. GENEBAL EQUATIONS 


In tke develcpmert of the equations which follow, itis 
assuned that these eguations can be validly used tc calcu- - 
late the flow over tke cone in either air or water. The 
equaticns were primarily developed by Shapiro [Ref. 8° for 
flow in air. = 
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II. FUNDAMENTAL EQUATIONS 


aw EACKCRCOUND 


Ike first ster necessary to describe superscric flow 
ever a cone is to calculate the thermodynamic freperties 
across a shcck wave. The basic research into the chance of 
thermcdyragic properties acrcess a shock wave in water was 
extensively concucted and reported upon in Underwater 
Explosicr Research [Bef. 1] during and just after Werld sar 
Il. Ike kest summarization cf tnese works can ke fcurd in 
Cole [Fef. 2]. 

The ~rigary source used as reference for the calculation 
cf tke Eydrcdynanic ysroperties of sea water at the ircrt of 
a shcck wave is the werk of Kichardson, Arons, and Halverson 
(Ref. 3}. They utilized graphical technigues, which were 
rTather crude and tedicus, te calculate the tnermrcdyranic 
data needed to describe the conditions of the sea water. 
Fuhs [Ref. 4] used tke work cf Richardson, et_al., [Fef. 3], 
to develop a ccmputer program for the HP4ICV hand-held 
calculatcr which efficiently calculated the same thermody- 
Tamic [froperties. Fuhs' [Ref. 4} programs provided the 
tasis for the FORTRAN subroutines, used in this work, waich 
calculate tke same tktermodynamic properties such as fres- 
sure, temperature, and density. 

Tke funcamental ecuations used to calculate the therzro- 
dynagic charges which cccur across an obligue shock frert in 
air have keen well krcwn fcr decades. These equations are 
excerticrally well described in Kinney and Graham [Ref. 5] 
and fcrm tke basis for the initial calculations fcr the 


superscnic flow fast a cone in air. 
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coe~dizersicnai gomertua eyguaticn states: 


ee ge 225 
P, + 0,V,/2 = py + pV, /2 (ee25) 


Now ncte tkat typical] values for pe, and vy for the precklen 
under ccrsideration are as fcllcws: 
P, = 1000 kgyn and Vo= 1500 m/s 


9 Fascals or about 11.1 kilo- 


Therefcre, om v 72 = 1.12 X 16 
kars. Typical values for Pi the upstream pressure, are on 
tne crder cf 1 tc perkaps a few tens of bars. Thus, it can 
ke safely assumed that F, can be ignored ccmpared tc the 
tpstreagz dyramic pressure. 

Eaving specified this jressure, the program of thkis 
thesis utilizes the subroutines developed by Fuhs [Ref. 4] 
to calculate the density and , Velocity components at a feint 
just ¢ccwrstream cf tke shock frent. Since the gecmetry of 
the velccity flew across tke shock front is the sare in 
€ither air cr water (see figure 2.4), the velocity ccmfo- 
rents detergzined in the first step can be used to "jump pdack 
across" the shock frert to the upstream side where tie free- 
streaa Mach rumoker a can te determined. 

After determining the conditions on the downstream side 
cf tke stock front, the program iterates along the strean- 
dines tc determine tke cone séemi-vertex angle. This itera- 
tion is ferformed using the Busemann graphical construction 
as for tte air case. The majcr difference between the water 
and air cases, after the initial shocked conditicns have 
reen deterained, is that the local speed of scurd cn a 
streaglire must be determined iteratively, for the water 
case, ky using Fuhs‘ [Ref. 4] subroutines tc detergine the 
density and pressure at each foint. Knowing the pressure 


and dersity, the syeed of 


scurd 
lated trem: 


in the water can te calcu- 
B, ,c,n ag 
fe ES (2.26) 
pe 
woich is derived fren the modified Tait eguation usinc the 
defiriticne cf the speed of scurd as: 
2 _ 2p 
ee 2 Sp Is 


These eguations, 


ccikined with the equations develcped in 
Secticn — cf this chapter, constitute 


all the eguaticns 
needed tc ccmpletely solve the supersonic flow over a ccné. 
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| TABLE i 


ccerical shock cr a sheck front 


| ()  Signifies a condition downstream of 
2 cenical sheck or a sheck front 
| 


( ). Signifies ccnditions at cone surface 
(). Signifies a ccmponent in r-direction 


() Signifies a component in w-direction 


> 
| 
£ | Nomenciature a 
: t { ee 
c Speed of scurd ee 
Ms k Ratio of specific heats fc 7c) | pose 
x Mach nuaber 
Me Mclecular weight of air | 
E Pressure 
{ L Racius in spherical cocrdinates { “ 
F Radius cf curvature of the 
| hcdograrph streamline | 
1 Temperature 
| v Velocity | 
Vax Maximum velccity for adiabatic flow :, 
B Angle the streamline makes with the shock ve 
Flane {same as the angleo) in air calculaticns ae 
: | 8 Ccne seni-vertex angle i oe 
: 8 Flew directicn | ek 
0 Mass density mo 
" o Shceck argle | 2S 
4 J w Angle in spherical cccrdinates j a 
: A Universal gas constant | oo 
. ( , Signifies a ccnditicn upstream of | 5 
| 
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A. EFBCGEAE LOGIC wo 


As discussed in Ckapter 1, the computer progran devel- a 
cped fer this thesis sas written in the FORTE&AN programming 
Janguace. Structured pregramming fractices have Leen 
folicwed throughout in that the program is divideg ito 
Flocks (cr acdules) cf code, each of which is designed to 
Fericrm a single calculation sequence. Many oi the agcdules 
have keen -laced inline rather than being writter as seéepa- 
rate stkrcutines or ftnctions. This was done prisarily for 
€ase cf use and understanding. In addition, however, mecst Bey 
modules were placed inline because the farameter list rae 
exchanced between the module and the main program would aave 
keen excessively lcng ctherwise. Regardless of whether the 
code is included inline or as a separate module, the fregran 
waS writter in a manner which ensured that the "“side- 
nt effects" problem disctssed Ly MacLennan [Ref. 9] did rot 
” cccur. In addition, the use of FORTRAN'sS CCMMON ccnstrict 
= has keen studiously avoided to prevent the aliasing precbien — 
z discussed Ly MacLennan [Ref. 93. A module hierarchy chart, 
whick shows the matcr modules of the program and their vee 
interccrrections, is c¢iven as figure 3.1. 

The flew logic cf the main computer program is illus- 
trated ry figures A.1 through A.7 of Appendix A. These flcw- yt, 
charts shcw that the initial part of the program (uf tec line 
60) is used tc initialize certain key variables and to Reece 
gathér irput values from the frogram user. Note that, as oe 
stated in the objectives for this proyram, all user input is 
_ verified tc ensure ‘that the values entered fall within 
-. Frescrited limits. In crder to ensure that input values are 
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witnixc limits, the sukroutine CHKINP is aiunvoked and the 
input values are passed as ,;arameters. The flew lcgic ci 
subrceutire CHKINE, whicn is self-explanatory, is ircluded as 
figures A.8 through A.10 of Appendix A. 

After receiving and verifying the user's ingut, the a 
Frogram tesgins the calculaticns reyuired to determine the aoe 
cone séni-vertex ancle. Frem line 60 to line 40C, the 
Froygram calculates the initial thermodynamic projyerties of. 
the water cr air, as the case may be, and calculates the 
changes which occur in the fluid properties as a result of 
the fassage of the shock front. As mentioned previcusly, 
these calculations fcr the air case are ratner straightficr- 
ward and, kecause of this simplicity, the computer ccde for 
the air calculation tas been slaced inline. For the water 
case, tke subroutines WSHOCK and EPV are needed to calculate 
the required thermodynamic properties. Subroutine WSHCCK 
calculates the thericdynamic freperties of the water at 
Foint 1 upstream of the shock front and at point Zz, jrst 
across tke shock frert on the downstream or shocked side. 
Appendix A, figures 2.12 and A.13, contains the flcwchart 
which describes the lcgic of this subroutine. Subroutine EFPV 
was ccried, with the author's rermission, from Fuhs [Ref. 4] 
and is fully describkec in that work. However, the FCRIRAN 


transiaticn of EPV is included in the prograr iistirg 
contairnreé ir Appendix C. 


Eavirg calculated tae initial thermodynamic fluid ccndi- 
tions cr tcth sides cf the sheck front, the program begins 
the iteraticn precess reguired to determine the cone seni- ee 
vertex argle. These calculations peyin at line 400 cf the a 


computer gregram. Tkis iteératica process makes use cf the 
Eusegann graphical integration technique discussed by 
Sharire [kef. 8]. First, the radius of the Busemarn curve 
for a feint J is calculated using the velocity at that 


Foint, tke streatline angle at that point, and the value cf 
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the argle cuega at that ;oirt (see figure 3.2). Next, tne 
center of the circle for tne Eusemann Curve at the point J 
is calculated using the velocity components at feirt J 
(whick are desigrated U and V), the radius calcuiated ir tae 
Frevicus step, and tke value of the angle omeja. Tais circle 
center ysrowides the scint from which the prcogram "draws" the 
arc used tc calculate the next foint on the cusemann curve. 
Finally, tke velocity and the Busemann velocity ccmgcrents, 
at the next foint, Jt#1, are calculated by "“swinying an arc" 
from the circle center calculated in the previous step. in 
additicr, the streazline angle at the next feint is 

calculated. 

The ;rogram next tests tc see if the cone surface has 
keen reacked (line 475). This test is conducted ry deter- 
ginin¢e if the aksolrte value cf the difference hetween the 
streamlire angle and the angle omeya is iess than a sfeci- 
fied test value (which is set to 1X 107%. If this difizer- 
ence is less than the test value, the cone surface has ceen 
reached. In this case, the frogram then calculates tke ther- 
nodynaric fluid conditions at the cone surface and the ccne 
semi-vertex angle ard dispflays the final results cf the 
Frojram. If the difference Letween the streamline angle and 
the angle cmega is creater than the test value, the ccne 
surface Eas not beer reached. ' In this case, the pregran 
calculates the therscdynamic fluid properties at the next 
point, J#1, then loops rack to iine 400 to kesin the 
Eusegann graphical irtegraticn f-rocess for the new feint. Ina 
the case of a water run, the subroutine WAIVEL is irvoxed 
just Ericr to lcoping back tc line 400. This sukrcutine 
calculates the thermcdynamic properties of the water at any 
given ;cint. Apfrendis A, figures A.14 and A.15, gives a 
flowchart which demcrstrates the logic flow witain this 
subrettine. 
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The ¢rcegram incluces a feature which allows the user to 


take repeated executicns of toe program without the reguire- 
Bent cf "reloading" tke frogrags kéetween each executicn. This 
feature was included fy soliciting a response from the user 
at line 16C€0 as tc whether heyshe wishes to make anctker 
executicr cf the program. An affirmative response causes the 
Frogray tc loop tack to statement 1 at the beginning cf the 
Initialization secticn of the main jprogram. A negative 
respcrse causes frogram termination. 

Nete that certair separately compiled functions included 
with this program have not teen flowcharted because their 
logic is sc straightforward and because the subprocrags in 
guesticn usually consist of cnly one or two lines of execu- 
table ccde. The furctions fitting into this categcry are: 
[TOK, BICL, and SOSAIF. 


Ee USEE INSTRUCTIONS 


As disctssed in the objectives listed in chapter 1, cne 
goal cf this thesis sas to ensure that the program presented 
waS e€asy tc use and nsaintain. Ease of use has been facili-~ 
tated ky including code which verifies that the input 
provided ky the user is within specified ranges. For 
example, in the sater calculations, the pressure at feint 2 
can Le ne greater tkan 100.0 kilcbars. This is due tc the 
fact tkat Fuhs' (Ref. 4) sukroutines, which are used to 
calculate tke thermcdynamic properties of the water, are 
kased cn the work of Fichardson, et. al., [Ref. 3], which 
cniy gives results up to fyressures of 100.0 kilctars. 
Therefore, while tke subrcutines could comgute the water 
conditicrs at pressures greater than 100.0 kilobars, it is 
uncertain wkether the values so calculated would te entirely 
correct. Fcr this reascn, the ranye of the input values has 
keen restricted. Ease of use is furtner facilitated by 


as 
dar 


Frovidincg clear, méaningfei cutput. Aprfendix B ccrtains 
Sample cttpets generated oy the ;,Lrogram for the Sumzary and 
complete print cpticrs of the program fcr roth water ana > 
air. lat 


Cne cf the sajor criticisms leveled by computer scien- 


tists against the FCFIRAN lasjuage is for its lack cr a 
reguirezrent for formally defining all variables used within " 
a precrag. Similarly, FORTRAN has an implicit declaration 
Folicy stereby any variabie whose first letter is I, J, K, 
I, & or N is implicitly declared to be of tyre Inteyer. In 
this fEregran, these two unfortunate characteristics of ’ 
FORTRAN are avoided ty reyuiring explicit declaration of all 
Frograz variables. 

Ease of maintenarce and ease of understanding cf the 
Frogran kave been ackieved througn the use of a literxally aes 
commented frograz ard through the use of relatively sean- 
dagful variable names. In this regard, the FCRIRAN lancuage 
is less thar desirable since it limits variable names tc six 
characters. Appendix C contains a fully documented listing 


wee 


cf the ccmwputer program ard all subroutines or functicns 
used ky the main procgram (cther than standard lirrary turc- 
tions such as sir). Tables VI, VII, VIII and IX, lecated in 
Appenéix A, contain lists cf all variatle names used in tne phe 
Hain ~rocran, their weaning, and their MKS (seter-kilccran- 
second) units, if appropriate. Table X contains a sinilar 
variakle list for the strroutine WSHOCK and fTakle (XI 
contains a variatle list fcr the subroutine WATVEL. ‘These _ 
two takles ére also lIccated in Appendix A. 

Finally, tais prcgram was developed usinj the FORTRAN IV 
danguace supplied by the ITEM Ccrporation as part cf their 
IBM 27Cyi023AP comptter syster which is the main ccmguter ‘1 
systeéea available at the Naval Fostyraduate School. As far as 
is knewr, no implenentaticn specific features have keen 


included in this pregran. Therefore, the program stculd oe 
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execute cn any systen which has a FORTRAN compiler. In créer 
to execute the frograa on the IBM 370 ccmputer syste avail- 
able at tae Naval Postgraduate School, the tolicwirg sters 
must ke acccmplished ir the crder givea: 


(1) The pregranz must te compiled using the comranc: 
FCRIHX CCNEFLOW 


This ccugard invokes the FORTRAN H Extended compiler, which 
is ar cptinizing, freduction run compiler supplied as part 
cif tke IBM comguter system. The proyram could alsc te 
compiled using another FORTEAN compiler such as the FOFRTEGI 
compiler. Note that the abcve assumes the frograa suyylied 
Ey this thesis has teen entered into a file wnich has the 
filename CCNEFLOW and whick is cf filetype FORTRAN. Ncte 
further that this step need only be performed cre tine 
Frovided errors do nct occur. 


(Z) Next, the libraries of standard subroutines and 
functicrs stpplied by the ccugputer center must tLe attacted 
to tke file. This is accomplished by issuing the fcllcwing 
coppard: 


GIOBAL TXTILIE FCRTMOD2 ACD2ZEEH 


Note that this step need cnly be performed once [rcevided 


abncexrmal terminaticnus dc not cccur. 


(3) In the FORTEAN language, a fcrmatted inputycutyut 
(1/0) statement regtires the definition cf an Iy0O device. 
For e€xargle, at most installations, by default, IyC device 5 
is used for input ard IyvyO device 6 is used for cutyut. To 
utilize the program ci this thesis properly, the fcllcwing 
commards mtst be issued to define the IvO devices tc te 
used: 
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FILELEF 06 LISK CONEFLOW CUTPUT 
FILELEF O07 TEnY (PERM) 


As a result of the akcve cciwands, all output frer tke 
Frogram till be written intc a file on the user's disk raed 
CONEFICK with a filetype cf CUIPUT. It is not strictly 
necessary tc issue tre FILEDEF 06 command shown akcve, cut, 
af nct, all outfut will be written at the terminal and will 
not ke available for printing. The second file defirition 
commard (FILEDEF 07) is absclutely reguired or the pregran 
will rect cperate. Ikis file definition tells the cczputer 
that IyG unit 07 is tke cowuputer terminal. In the p;regraa 
ci this thesis, all input from the user is reguested fren 
1/O urit Ci. Therefcre, if unit 07 is not defined, tae 


Frogram cannot receive any input. 


(4) Finally, the program can be executed Ly issuing 


the ccruand: 
ICAL CONEFLOW (START 


when frocram executicn is completed, a cryptic messace of 
the fer R; T=0.01,0.01 16:45:00 will be displayed at the 
terminal. This is siajly a message from the ccaputer's cper- 
ating system indicating that the task just requested sas 
teen cctypleted. The user can now utilize the system's cper- 
ating ccugands to review the output from the executicn of 


the yicgiaz. 


If the user desires to execute tae frogram again, return 
to step (3) of the akove procedure. If the cutput frea the 
Frevicus executican is ne longer needed, the same FILEDEF 06 
ccumand can te issued. However, if the previous execution's 
cutrpeut is reeded, tke filetype of the FILEDEF 06 ccrmand 
should ke chanyed (e.g. change OUTPUT in the above ccmrmand 
to CUIELIZ). 


IV. PRCCRAM CALCULATION RESULTS 
A. FERCGBAM RESULTS FCR THE CALCULATION OF CCNICAL FICH IN 
AIK 


Ike pregram accurately calculates the flow over a ccne 
in air. Numerous prcgram executions were made for tke air 
case at various upstream Mach numbers and for varicus shcck 
angles. The results crtained from these program executicns 
are sumgarized in Takles II and III. The numbers reduced Ly 
the Fregrag were ccopared for tne variables showr int the 
tables tc the tables given in Kinney and Graham [Fef. §], 
against the yrafhs civen by Shapiro [Ref. 8], ard tc the 
tables produced by kcpal [ Bef. 10]. The vaiues given in 
thesé scurces are ircluded ir Tables II and III, woere 
apprcyriate, alcng with the calcuiated results. As can be 
seen frcxr taese two tables, the calculations made Ly the 
frograr cive yuite accurate results (less than a 1% exrcr in 
nost case€s). Certairly, the program is much more accurate 
than cne's ability tc read the graphs jyresented in Shapiro 
[Ref. 5. 

Eased cn these ccafariscns, it is believed that the part 
cf the frogram which calculates supersonic tlow over a cene 
in air is accurate. Mhe importance of this fact is that, by 
feelirg confident that tne procedure followed for the air 
case is valid ard that this frocedure nas been ccrrectly 
juplenented in the froyramming language, it is safe to 
assume tkat the Eusegann calculaticn procedure utilized ter 
the air case can ke accurately applied to the water case 
Frovided the water ccnditic:s at each point are calculated 
correctly. 
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As discussed if the Introduction, this thesis hac 
certair gcals wnick it was desired to achieve. It is 
relieved that these coals have all been successfully met. 
The gregraz of this tkesis has Leen successfully translated 
from Easic into FORIBAN as desired. In the translaticn of 
the criginal prceyran, modern structured programmirg frac- 
tices nave been follcwed tc the greatest extent pessible. 
Finally, tke ¢regraga freserted is guite "user-friendly" and 
is very well-documented. 

As was demonstrated in the results chapter, this gfregran 
Calculates correct results fcr each cf the air cr water 
cases. It is mentioned here that, as pointed out freguently 
Ey ccuputer scientists, it is virtually impossible tc test a 
Erogram fcr all possible cases. Therefore, no fregram is 
entirely error-free. It is Lelieved that the prograr cf tkis 
thesis is as free of errors as is possible without exhaus- 
tive, tine-consuming and extremely expensive testing. 

As a result of having get the three goais mentioned in 
-the first paragraph, it is kelieved that the prograr cf this 
thesis will frovide an excellent werking tcci tor 
researckers in the field of shaped-charye jet penetraticn in 
water. Tbe fact that the precgram is well-documented wili 
jake any ucdificaticrse or extensions to the program, shovwld 
that ke reguired or desired, easy to perform. Further, sirce 
the frcegram is e€asy to use, cnly a cursory kncwledge of 
computers i¢ needed in order to utilize the frogran. tMhkese 
features are regquiregents fcr any computer program which is 
to be used Ly scientists as_a_tool. Too often, computers and 
their gpregrams reyguire that the people who simply desire to 
use tkeir capabilities must learn a great deal of detail 


in 
®W) 


39zROK UT MOTT OTUOSIadAS FOJ YasTOTJJa0j Heig Sa Jaqunn YOuN HerxIZS3AaTG Evh 2INGTT 


> 0@='1 YIGWAN HOW WYIYLSIING 
0°S S's O°b S°e o’e 


t 00°0 
S0°O 
co ae) 
St°O 
02°0 


w 


it 
HT 


ATU 


a 


Os.°0 
S£°0 
0b°O 
Sb 0 
0Ss°0 
ss°0 
03°0 
s3°0 
02°0 
$Z°0 
08°C 
S8°0 
06°0 
$6°0 
00°t 


Ht 


LN3T9144309 SYNC 


YLT, 


z 
ALT 


Yalawuddd SU FTONY SNOD HIM 
Y4EWNN HOBW WYdYLSA4Y4 SASYIA 
INSTOT343509 9vyC 4O LO Nd 


ee ee ana nar ee, 


51 


ya ee 
SE re Rel CURR ir rial roe ric SU ie ante (Ma cate ie OL died BO Se 


*1a}zeM UT MOTH DtuosIadns FOF WANN YOeW eoeyINS SA TVIINN Youy tearysaarg 


7° Danbry 


3,02 =1 YIGWAN HOUW WYIYLSIGYG 
Q’°sS G*b 0° GE ones S°z o°? S°l a°t 


(o> es 2 n° oO © 0, OEE 3 


YSEWNN HOW 


YATIWUNdd SU A TENU NOU HLIM 


MAGWAN HOBW WESYLSI4Y4 SASYSA 
YIEWAN HOWW SOYIYNS 3NOD 40 LO ld 


(re oe cr wr ere ee eae ne Yee Se NE ry EE TS SY SF EU et ee 


Oo oO oO OQ oOo Oo oO oO oO oO 2 
FOWsYNS 3NOI 


i rn rr i re ee a a nn i a ne ee = 


—— 


°79}zem UT AOTY OTaosraIns O¥ atTHuy yous SA TSQCHNN Youn TedIASIAIYT 


2 0¢=11 YAGWNN HOBW WeaYLSIIY4 
0°S : sf ove sz 


WILE SY SONG INOO HITM 
UIGHNN HOU WUTeLSTaas SASHIA 
TINY HOOHS JO LOTd 


b°h axrnbry 


rn na i a nn nn i a ee ee 


a°s 
oot 
o’st 
0°02 
0°S2 
0°OL 
O'St 
0° Ob 
O°Sb 
a°os 
O0°SS 
0°99 
0°S9 
0°04 
O°SZ 
0°08 
0°S8 
0°06 


AOOHS 


JTSNY 


ee ee ae cee ee Se ee ee ee ce me ee ce ee ce ee a te ee ee ee See a 


49 


computer system is being utilized at near its capacity, the 
same resrpcrse tine Car increase up to about Z to 3 Linutes. 
This result, however, is Symptomatic of any computer systen 
fteincg utilized at near capacity and iis not a result ora 
reflecticn cf the precram's structure or design. 


4E 


greater than 5.0 (a converience for ease of presertation 
cnoly). The pre-precessed data points were passed tc an 
interseiaticn and smoothing program develozsed fry Frarke 
[Ref. 12° which used a guadratic Shepard's method tc smocth 
the cata points generated. After smoothing, the data was 
fost-;rocessed to remove erroneous interpclaticn feirts 
which were generated because the original data was irreg- 
ular. Finally, the processed data was passed to a FORTRAN 
Frogrags developed ry the author of this thesis thich 
utilized tke Contcuring feature of the DISSPLA grarhing 
routines of the ISSCC company. The results of this extexusive 
Frocéssirg are the graphs presented in this thesis. FEecatse 
cf tke irregularity of the original data (which causes 
difficulties ia tke quadratic Shepard's interjclation 
gethcd), and because of the smoothing reguired Fy the 
contouring routines cf DISSFLA, the graphs presented are 
somewhat rough and skculd, therefore, only be used for "back 
c£ tke ervelope" calculaticns. Accurate results are 
Frovided Ey the CONEFLOW prcgram which should be utilized 
for mere precise work. 

Finally, one of the primary reasons for translating the 
crigiral Basic frograms into a higher-order language (in 
additicr tc the desire to make the programs more accessible) 
was tc speed the execution cf the proyran. In this regard, 
tne werk cf this tkesis has nore than accomplished this 
result. It aust be menticned here that the pregrat will 
execute extremely fast, considering all the iteraticns 
needed, prcevided the computer system in use is net heavily 
in use at tke tize. For example, at the Naval Postgraduate 
Schocl], when the ccrputer system is relatively free, the 
Frogram cgerates so ctickly that one execution requires less 
than £ seccnds (from the time cf the last user input tc the 
time cf the reguest by the fregram for an indicaticn of 
whether ancther execution is desired). Conversely, when the 
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This ccrclusion is true provided the upstream Mach nuaker is 


held ccnstant by varying the pressure oehind the shcck frcnot 
throuch the mechanisno cf the Eressure aultiplicaticn factcr. 
Having determined tlkat tae caiculations can be ccnéucted 
independently of the temperature, plots of the variakles of 
interest ir water sere made in a manner similar tc the 
craphs fcr air shown in Shagirc [Ref. 8]. 

Finally, having determined that the calculaticns ccula 
tke mace independently of the upstream temperature and tkat 
these calculations sculd be accurate (within 1% for ccne 
angle and 123% for pressure), it was decided to plct varicus 
Farameters cf irterest for the flow over a cone in water. 
These craphse are presented as fcllows: 


(1) Figure 4.1 is a plot cf the freestream Mach Numter 
versus the shock angle with the cone semi-vertex 
angle as an entry fparameter. This grafh is the 
water analog to Figure 17.7.{a) of Shapiro [Ref. &]. 


(z) Figure 4.2 is a plot cf the freestream Mach runter 
versus the Mach numker at the cone surface with the 
cone semi-vertex angle as an entry parameter. Ttis 
cragh is the water analog to Figure 17.7. (c) of 
Shazyiro [Ref. 8j. 


(3) Finally, Figtre 4.3 is a plot of the freestream Mach 
Lumrer versus tne drag coefficient (c ) with the 
ccne semi-vertex angle as an entry parameter. This 
cragh is the water analog to Figure 17.7. (f) of 
Shagiro [Ref. 8]. 


It shculd be nected here that the graphs develcred and 
Fresented in this tkesis were obtained through repeated 
executicns of the CONEFLOW frogram, which generated approxi- 
mately 3€00 data pceints per graph. These data points were 
Fre~;Erecessed to rencve entries with a Mach nugker of 
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Farareter, does not exist. Since water is net a thernrally 
and calcrically perfect fluid, it is not possitrile tc slot 
universal] craphs of tte variatles of interest. 

Ic gain insight into the sensitivity of the supersonic 
conical flcw in water to the upstream conditions, a study 
was ccenducted in whick calculations were performed tc deter- 
Mine swhetlLer the cone semi~vertex angle is relatively inde- 
fendert cf the upstream temperature or whether varicus 
Frograp executions are required to Snow the variaticn of 
cone se€ni-vertex angle with water temperature. 

Ike cGéeterminatior discussed above was made Ey hcelding 
the uystreaz pressure, the shock angle, and the urpstrean 
Mach nuagker constant. Then the upstream temperature sas 
varied tc determine the variation in the ccne semi-vertex 
angle. In crder to kcld M, constant, the fressure multifli- 
caticn factcr, described in Chapter II.D, was varied until 
the Mach number upstream for the new temperature matched 
(withir reasonable accuracy) the upstreau Mach cumster tor 
the crigina] temperature. Note that “reasonable accuracy" 
for aatchinsg of the upstream Mach nuagbers meant matching the 
pumcerse to the fourth deciazal place. Table V prevides the 
results freu these calculaticaos. As can be seen fren this 
table, the cone semi-vertex anyles calculated fcr the 
varicus temyeratures are all cf approximately the sare value 
(withir 1%). The fressure dcwnstream of the shcck wave 
varies Ly 13% as can fe seen in Table V. While there is small 
variaticr tetween tke values for the cone semi-vertex 
anyles, it wust ke remembered that, for practical furpeses, 
these variations are very small and can be safely ignoied in 
examining the flow cver a cone in water, Thus, it was 
concluded ttat the caiculaticn of the cone semi-vertex angle 
in water can be conducted independently of the water tenper- 
ature {at least for the accuracy reguired for the calcula- 
tion cf the movement cf the metal jet from a shaped-charge). 
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liezaganr and Rcshko [Ref. 11] demonstrate tkat tne 


thermal énd calorac eguations cf state are thermocyranzically 
related ty: 


me (4.6) 


Intrceducing equation 4.2 itte equation 4.6 results in the 
Light-hard side of eéguaticn 4.6 becoming “egual tc 2eL0. 
Thus, fcr a therzally frerfect gas, tne enthalpy is indeépen- 
dent cf the pressure, and, hence, enthalpy is a furction 
cnly cf temperature. ‘Therefore, a necessary conditicr for 
eqguaticn 4.5 to fe valid is to have a thermally perfect gas. 

Fer water, Richardson, et. al., [Ref. 3], demcnstrate 
that tke Leat capfacities are functions of the temperature of 
the water. Thus, ecuations 4.4 and 4.5 are not valid fer 
water, and water is ret a calorically perfect fluid. 

If a fluid is Ectk theroally and caiorically ;yerfect, 
the ccnditicns acrcss a norgal Shock front depend orly on k 
= Cp/c, and on the freestream Mach numker. Conversely, if 
the fluid is not thermally and calorically perfect, addi- 
tional variables must be specified in order to define the 
shock ccrditions. Extending the argument, one can state that 
the suyerscric flow cf a thermally and calorically ferfect 
fluid over acone is a furcticn only of the heat capacity 
ratic, k; the freestream Mach number, Mi and the shcck 
front angle, cg. As a result, a single graph is sufficiert to 
represent all flcws cver the cone. This is the conditicr fer 
a perfect gas. 

Ir ccntrast, when the fluid is not thermally and calori- 
cally perfect, the sclution fcr supersonic conical flew is 
Gependent upon variables cther than k, Mis and o . 
Consecuertly, a universal graph of the shock angle versus 
the freestream Mach rugber, with the cone angle as an entry 


£& 
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which describe the cone flcw under all corditiors as is 
shown ty Skapiro [Ref. 8]. whoetner universal results for 
tne water case could te oktained was not so cléar; this 
Feint is ncw discussed. 

Liegygsarn and Rostko [Ref. 11] discuss the general ther- 
mwodynazics of fluids and introduce the concert cf a 
"thergal" ecuaticn cf state. In yeneral, the thermal equa- 
tion cf state for any fluid is given by: 


£(EF,%¢1) = 0 (4.1) 


For a perfect gas, the thermal equation of state is given 


(4.2) 
where F =A /M,. Anotker equaticn which relates the therimccay- 
ramic variakles e, ve and Tis the "caloric" eguatica of 
state, which, in a general form, is given fry: 


£(e€,v,1) = 0 


A caicrically perfect gas is defined by: 


(4.5) 


where € is internal erergy and h is enthalpy. Thus, the heat 


capacity (either c, or c) is a constant for a calcrically 


perfect cas. 
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Eo EROGEAN RESULTS FCR THE CALCULATION OF CONICAL fFIC&# IN 
WATER 


As mentioned in the previcus secticn, the calculaticn 
os Frocecure was verified for the air case by compariscn wita 
known results. This ccmparison showed the method used ir tke a 
calculaticns was correct. ‘he next step taken was tc deter- - 
mine if the subroutines used to calculate the theracdyranic 
Froperties cf the water gave results ccmparable tc thcse of 
Richardscn, et. __al, [Ref. 3], and Fuhs [kef. 4]. 

Calculaticns were fpexrformeé which gave results which could -— 
Ke ccapared to Table JI of Fuhs [Ref. 4]. Tatle IV freserts . 

a Ccujyariscr betseen the numeric results calculated [Ly tke 

Frogram and those given by Fuhs. AS can be seen, these 
results are, within reasonable error, remarkably similar. 
This is net unexpected since the subroutines used [Ly the 

Wain y~rocraz to calculate these variarles are essertially 

direct language translations of Fuhs' frograms. Therefore, 
the results should be similar. a 

Havirg verified that the procedure utilized ir the ‘i 
calculaticn of the cone semi-vertex angle was ccrrect 
(thrcugh the air results) and haviny verified that tne 


FrOgram correctly calculated the thermodynamic precperties cf a 
the water at any peint, it was believed that the fregran fe 
could ke executed fcr the water case, for various initial 


conditicrs, with certainty that the results so calculated 
would ke accurate. Ecwever, a guestion arose as tc whetter 
the calctlations perfcrmed in the water case were indepen- 2 
dent cf the upstream thermodynagic properties of the water. 
In ait, tke only guantities needed to calculate the ccne 
semi-vertex angle are the upstream Mach number (4)) and the 
shock angle (6 ) (€.y- see figure 17.7(a) of Shapiro ". 
[Ref. 8] ). Thus, for air, the calcuiaticn of the ccne 
semi-vertex angle is independent of the upstream temperature ~ 


cr presstre, and, therefore, universal curves can be drawn 
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about tke computer in order to use the tool. In téis 
Frogram, ar understanding of FORTRAN is nect reguired in 
crder tc utilize the jrogram cr its design. Unfortunately, 
due tc tke cperating system cf the IBM 370 computer syster, 
the same carnot Le said about the steps reguired in crder to 
actually use the pregraz. Scme familiarity with the ccryputer 
systea in use at the users location will re a necessity in 
crder tc cperate the [regrag ccrrectiy. 

Finally, as mentioned in the Introduction tc tkis 
thesis, tke program till provide an excellent test case and 
ccomfariscn nodel for a computer program which mcdels the 
actual ficw of the metal jet from an explosive sharped-charge 
fired thrctgh the water. The actual situation is a Elurt- 
nosed, sxatter tnaan a sharp-pcinted conical, flow precblen 
which is ruch more difficult toc solve. Therefore, a kncwn 
soluticn and methodolcgy of sclution for the easier precblen 
is a necessary first step to the solution of the larcer 
Froblem. It is believed that the program of this thesis will 
serve as this necessary first step. Further, it is kelieved 
that, wken the equaticns solving the actual flow prokles are 
develcred in their final fora, the proyram of this thesis, 
due tc its ease cf acdification, will serve as the pregras- 
fing mcdel for the pregram which calculates the biurt-nosed 
flow freckled. 


EROGR 


This appendix cortains the logic flowcharts of the main 
Frogragz and its subrcutines (i.e. those which have ret been 
descrited elsewhere cr waich are not part of standard 
computer center libraries). The language is kept rcatter 
general so the cverall logic cf the proyram can te demcn- 
strated. If the reader désires to know how a particular 
logic sequence is implemerted, he need only refer Girectly 
to tke fregram secuent which the logical flowchart is 
descritirg. The flicscharts included in this ajgpendix are 
listed Lkelcw: 


(1) Ihe Main Procram Flowchart consists of Figures 4a.1 


thrcugh A.7. 


(z) Function CHKINP consists of Figures A.8 
thrcugh A.10. 


(2) Sukroutire DEFANG ccnsists of Figure A.11. 


(4) Sukroutine WSEOCK ccnsists of Figures A.12 and A. 13. 


(5) Sukroutine WAIVEL ccnsists of Figures A.14 ard A.i15. 
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APFENIIX B . 
SAMPLE FRINIOUTS nee 
This apyzendix presents copies of the output From varicus ee 
runs of the computer fprograz listed in Apfendix Cc. The — 
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(2) A Summary Print of a Water Run 
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